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Abstract. Application of the standard canonical quantization rules of quantum 
field theory to macroscopic electromagnetism has encountered obstacles due to 
material dispersion and absorption. This has led to a phenomenological approach to 
macroscopic quantum electrodynamics where no canonical formulation is attempted. 
In this paper macroscopic electromagnetism is canonically quantized. The results apply 
to any linear, inhomogeneous, magnetodielectric medium with dielectric functions 
that obey the Kramers-Kronig relations. The prescriptions of the phenomenological 
approach are derived from the canonical theory. 

PACS numbers: 42.50.Nn, 12.20.-m 
1. Introduction 

Macroscopic electromagnetism is an indispensable part of theoretical physics, providing 
an accurate description of a large variety of light-matter interactions [TJ [2] . The ability 
to encompass the electric and magnetic responses of media within macroscopic dielectric 
functions, without reference to the microscopic material structure, is of great practical 
importance. A microscopically complete treatment of the interaction of light with most 
real materials is currently impossible, and even if the computer resources necessary for 
such a task were available, simple analytical calculations using the macroscopic theory 
would still be adequate in many circumstances of interest. 

But the source of the strength and versatility of macroscopic electromagnetism, 
namely the removal of microscopic material degrees of freedom, has been seen as a 
disadvantage when it comes to quantizing the theory [31 H]. A rigorous quantization 
procedure would require a Lagrangian and Hamiltonian formulation of the theory, 
followed by the standard canonical quantization rules. Such a canonical formulation of 
macroscopic electromagnetism has so far been thwarted by the presence in all media 
of dispersion and absorption, an obstacle that cannot be avoided if one considers 
all frequencies of electromagnetic waves for which the macroscopic theory is valid. 
This difficulty has led to a phenomenological approach to macroscopic quantum 
electrodynamics (QED), in which no canonical formulation is attempted (some details 
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of this scheme will be described below; for a comprehensive description see [2J E] 
and references therein). The lack of a canonical basis for macroscopic QED is 
ameliorated by the results of calculations that introduce a simple microscopic model for 
a dielectric [U [7J [8j [9] . Huttner and Barnett [I] showed that the canonical quantization of 
the resulting coupled light-matter system leads to results that can be expressed in terms 
of a homogeneous electric permittivity, derived from the matter degrees of freedom, 
at which point the model reproduces relations introduced in macroscopic QED on a 
phenomenological basis j5j E]- The simple microscopic model has been extended to the 
case of an inhomogeneous dielectric [7J [8] and to a magnetodielectric medium [9], again 
producing results at the level of the dielectric functions that agree with the prescriptions 
of the phenomenological theory [5JE]. An essential ingredient of the microscopic models 
is the introduction of a reservoir in addition to the electromagnetic and matter degrees 
of freedom, in order to incorporate the dissipation that is necessarily present due to the 
Kramers-Kronig relations. 

The drawback of the microscopic models as a justification of the phenomenological 
approach is, of course, their simplicity. As a classical theory, macroscopic 
electromagnetism is experimentally proven to be valid in the case of highly complicated 
materials, and a rigorous quantization would ideally also proceed on a macroscopic 
level. This requires a canonical formulation and quantization of macroscopic 
electromagnetism, a task that is performed here. The resulting theory is valid for all 
linear, inhomogeneous, magnetodielectric media with dielectric functions that obey the 
Kramers-Kronig relations. The treatment is macroscopic, in that the electromagnetic 
fields are shown to obey the macroscopic Maxwell equations for an arbitrary linear 
medium, but a reservoir is necessary to account for dissipation, as in the microscopic 
models. Indeed, it is clear from the outset that a reservoir must be included in addition 
to the electromagnetic fields, since otherwise one would obtain a Hamiltonian involving 
only the latter, and thereby a conserved electromagnetic energy in dissipative media. 
But the necessity of including a reservoir does not imply that additional dynamical 
degrees of freedom for the medium must also be introduced; only the electromagnetic 
fields and the reservoir are required. This last point has been recognized in some previous 
treatments of QED in media [101 El [12j [13j [Hi [15]; we here give a comprehensive 
derivation of general, linear macroscopic QED from electromagnetic fields coupled to 
a reservoir. Due to the common feature of a reservoir, the details of the canonical 
macroscopic theory have much in common with the microscopic models [H [7J EJ E| . In 
fact, it is no more difficult to canonically quantize macroscopic electromagnetism than 
to quantize the simple microscopic models. But the canonical theory is a greater reward 
for a similar amount of labour, as it is applicable to any medium where the classical 
macroscopic theory is valid. 

The postulates of the phenomenological theory [HI E] that serve as its starting 
point are here derived from canonical quantization. We thereby remove the need for a 
phenomenological approach since the canonical theory has the same range of validity 
(in contrast to the microscopic models). As well as providing a foundation for the 



Canonical quantization of macroscopic electromagnetism 



3 



wide range of important results that have been obtained from the phenomenological 
prescriptions [3 [6], the canonical theory will allow investigation of issues that require 
a Hamiltonian and the related canonical machinery to be properly addressed. An 
important example where the benefits of a canonical theory are desirable is the Casimir- 
Lifshitz effect [161 El HI HE] , the phenomenon of forces on bodies due to electromagnetic 
zero-point and thermal fields. The standard Lifshitz theory [T71 [HI [19] of this effect, 
being macroscopic and lacking a canonical foundation, suffers from the drawbacks of 
the phenomenological approach to macroscopic QED in general. This has led to doubts 
about the status of Lifshitz theory as a proper quantum treatment; the lack of an 
intelligible Hamiltonian has been criticizedjj and it has recently been argued that the 
theory is essentially classical [20]. The canonical macroscopic QED developed here can 
provide a rigorous foundation for Lifshitz theory. 

The action principle from which all the results follow is written in Section [2] and 
shown to lead to the macroscopic Maxwell equations. In Section [3] the standard rules of 
canonical quantization are applied and the Hamiltonian is derived. The Hamiltonian is 
diagonalized in Section H] and the electromagnetic field operators are written in terms 
of the diagonalizing operators; at this point the postulates of the phenomenological 
approach emerge. 

2. The action and the field equations 

In the canonical formulation of electromagnetism [21] the dynamical fields must be taken 
to be the scalar potential <f> and vector potential A, related to the electric and magnetic 
fields by 

E = -V0-<9iA, B = VxA. (1) 

The relationship between fields in the time and frequency domains is, for the example 
of the electric field, 

1 f°° 

E(r,t) = — / dw [E(r, w) exp(-iwt) + c.c.1 . (2) 
27r Jo 

In limited frequency ranges where losses are negligible the action for dispersive 
macroscopic electromagnetism can be written as a functional of and A only [22J. 
Our concern here, however, is with the theory in its full generality, where all frequencies 
are included and the dielectric functions obey the Kramers-Kronig relations [21 1]- This 
means that in addition to </> and A the Lagrangian must include dynamical variables that 
account for absorption and dissipation of the electromagnetic energy. We incorporate 
dissipation in exactly the manner used in the microscopic models [UCOIE], by introducing 
a continuum of reservoir oscillator fields. As we allow for a magnetic as well as an 
electric response in the macroscopic medium, the reservoir will be a set of two such 
continua of oscillator fields, which we denote by X aJ (r,t) and Y u (r, t), where u > 
is a continuous label giving the frequency of oscillation. The dynamical variables are 



| I am indebted to Gabriel Barton for properly stressing this deficiency of Lifshitz theory. 
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thus 0, A, X u and Y w , and we must write an action that is a functional only of 
these fields and whose dynamical equations are the macroscopic Maxwell equations 
in a general inhomogeneous magnetodielectric medium. As there are no dynamical 
variables representing the medium, it must be included only through the coupling of the 
electromagnetic fields with the reservoir. 
The required action is 

S[<f>, A, X,, YJ = S cm [<p, A] + S X [XJ + Sy\Y u ] + S int [<p, A, X w , Y w ], (3) 

where S em is the free electromagnetic action: 

S em [0, A] = y y d A x Q^E ■ E - B ■ b) , Kq = 1/fio, (4) 

Sx and Sy are the actions for the free reservoir oscillators: 

S X [XJ = \Jd 4 x J™ da; (^X, ■ d t X w - oj 2 X„ • Xj , (5) 

S Y [YJ = ~ y d 4 x y°° da; ($Y W • <9 4 Y W - a; 2 Y w ■ Yj , (6) 
and Si n t is the interaction part of the action, coupling the electromagnetic fields to the 



reservoir: 



//*oo 
d A x J du [a(r, o;)X w • E + 0(r, a;)Y w • B] , (7) 

2k 



7T 



1/2 



(8) 



-,1/2 

7T 

In the coupling functions ([8]), £i(r, a;) is the imaginary part of e(r,co), the relative 
permittivity of the desired macroscopic medium, and Kj(r,u) is the imaginary part 
of k(t,uj) = /i(r,a;) _1 , where fi(r,u) is the relative permeability. For dissipative media, 
Si(r,u) > 0, whereas Ki(r, co) < because //i(r, co) = — Ki(r, uo)/\k(t, oj)\ 2 > 0. The 
coupling functions jS]) are thus real and so is the action ([3]); this will ensure that the 
Hamiltonian of the quantized theory is hermitian. The entire dielectric functions e(r, co) 
and k(t,u) of the medium, not just the imaginary parts, are specified by (jSJ) because 
the real parts are given by the following Kramers- Kronig relation |U [2] : 

£R,(r, co') — 1 = — P f du — ^ ' } , and similarly for k(t, co). (9) 

7T J Q OJ 2 - to' 2 

It should also be remembered in what follows that £i(r, co) and Ki(r,co) are odd 
functions of co, whereas e^(r,co) and kh(t,co) are even, so that e*(r, co) = e(r, — co) and 
k*(t,u) = k(t, —co), which is a consequence of the electric and magnetic susceptibilities 
being real pp. We have assumed the medium is isotropic, with scalar dielectric functions 
e(r,co) and k(t,co); anisotropy can be included by obvious modifications, but this is not 
done here as it would further tax an already burdened notation in what follows. 

The coupling of electromagnetic fields to a reservoir has been considered previously 
in PH m US H31 H US]- Here we P rove that the specific action ©-© gives the 
macroscopic Maxwell equations for the linear, isotropic, inhomogeneous medium with 
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dielectric functions e(r,uj) and k(t,uj). Variation of 
gives, respectively, 



), A, X w and Y w in the action 



/■OO 

e V-E+ / dwV-[a(r,w)X u ] = 0, (10) 
-k V x B + e d t E + / dw {a(r, cu)<9 t X w + V x [/3(r, w)Yj} = 0, (11) 



-d 2 X w - w 2 X w + a(r, w)E = 0, (12) 
-%Y U - w 2 Y w + /3(r, w)B = 0. (13) 

In order to find the independent equations satisfied by the electromagnetic fields we 
must solve ( fl2|) and ([13]) for the reservoir fields X w and Y w in terms of E and B, and 
substitute the results into ([10]) and ( ITT]) . To solve (TI2]) we write it in the frequency 
domain using fl2]) and obtain 

w ,2 X aJ (r, uj') - w 2 X w (r, u/) + a(r, w)E(r, uj') = 0. 

The general solution of (Hi 
a(r,w) 



(14) 



is 



X 



r, uj 



r E(r, u/) + % - uj')Z w (v) + 6(uj + w')Z* (r) 



(15) 



- (u/ + iw0 H , 

where Z^(r) is an arbitrary function independent of uj' that gives the solution of the 
homogeneous equation (( [12]) with E = 0), and + is an infinitesimal positive number 
introduced to deal with the pole at uj' = uj that arises from dividing ([14]) by uj 12 — uj 2 . The 
prescription chosen in ( [15]) for dealing with the pole amounts to a choice of boundary 
condition for the reservoir field X^; in the microscopic models [U E], which also 
include a reservoir, similar poles are encountered in solving the coupled system of field 
equations and a choice is required that leads to the desired class of solutions. Note that 
the prescription for the pole, and the form of the solution of the homogeneous equation, 
must respect the frequency-domain property X* (r, u/) = X CJ (r, —uj') that follows from 
the reality of the field X w (r, t) in the time domain; this property holds in (ITS"]) . In 
transforming (JTSj) to the time domain the integration over frequency is determined in 
terms of principle values [23] by 



uj 2 - (uj' + iuj0+) 2 
1 

~~ 2uj 



1 

2uj 
1 



1 



UJ — UJ' 



+ 



1 



ky0+ uj + uj' 



io/0+ 
1 



UJ 



UJ 



- + sgn(uj')i7r5(uj — uj') + P + sgn(a/)i7r 5 (uj + uj') 



17T 



r2 +—[5(uj-uj')-5(uj + uj')] 



UJ + UJ 

(uj > 0). 



(16) 



ur — uj' 

The restriction to uj > in the last line is appropriate because uj is positive in ([15]) 
(see the action (J3j) — (EJ) ) - Note again that after insertion of (TIB"]) in ( [15]) the property 
X*(r, u/) = X w (r, — uj') is seen to hold. Using (jSJ) in the case of X w , we obtain from 

am and am 



i r 

2tt Jo 



duj' 



air, uj 



UJ + 



UJ 



^E(r,a/) 



exp( 



-ia/t) 



c.c. 
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+ 



ia(r, oj) 
4u 



E(r, to) exp(— hut) + c.c. 



+ 



2n 



Z w (r) exp(— iujt) + c.c. 



(17) 



To evaluate the X w -dependent terms in the field equations I jTOj) and (Uljl we require the 
following integral, which we evaluate using (JS} and the Kramers- Kronig relation fl9]): 



do;a(r,o;)X w (r,t) = — P 



7T- 2 



do; 



da;' 



u 2 — a/ 2 



E(r, a;') exp(— ia/t) + c.c. 



£o f°° 

H / da; [iei(r, a>)E(r, u) exp(— iut) + c.c] 

2tt Jo 

1 f 00 

+ — / da; [a(r,a;)Z a; (r) exp(— ia;t) +c.c] 
27r Jo 

£o /"°° 

— £oE(r, £) H / da; [e(r, a>)E(r, a;) exp(— iut) + c.c] 

27r Jo 

1 /"°° 

+ — / da; [a(r, a;)Z w (r) exp(-ia;t) + c.c] . 
2tt Jo 



The second term on the right-hand side of the final equality in (fT8l) is the electric 
displacement D(r, t): 



D(r,i) = — - / do; [e(r, a>)E(r, a;) exp(— lut) + c.c. 

Z7T 



(19) 



Equation ( {TBI for Y w can be solved in an identical manner to that used to obtain ( jTTl) ; 
denoting the arbitrary function, corresponding to Z w (r) in (fTTl) . by W w (r), the solution 
for Y w (r,t) is 



Y w (r,f) = r&u' [l^^B(r,a/)exp(-ia/i) + c.c. 
2,71 J \_u z — UJ U 



+ 



B(r, uj) exp(— iut) + c.c. 



— W w (r) exp(-iwt) + c.c. 

Z7T 



• (20) 



To evaluate the Y w -dependent terms in the field equations (jTOl) and (jTTjl we require the 
following integral, which is analogous to (ITS]) and is evaluated in an identical manner: 

f°° K o /"°° 

/ do> /3(r, o;)Y w (r, t) = K B(r, t) / do; [re(r, o;)B(r, a;) exp(— iut) + c.c] 

Jo 2tt Jo 

1 r°° 

+ — / da;[/3(r,a;)W w (r)exp(-ia;t)+cc]. (21) 
2vr Jo 

The second term on the right-hand side of (1211) is the negative of the magnetic H(r,t) 
field: 



H(r,t) 



2tt 



do; o;)B(r, lo) exp(— io;t) + c.c. 



(22) 



Substitution of f|T8l) . (jT9l) . (pTj) and (1221) into the electromagnetic field equations (jTOl) 
and (TTTj) yields the macroscopic Maxwell equations 

V • D = a, (23) 

VxH-^D=j, (24) 



(25) 



(26) 
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where the free charge density a and free current density j are given by 

1 f°° 

cr(r, t) = V- / du [a(r, u;)Z w (r) exp(— iut) + c.c] , 

2f io 

1 f°° 

j(r,i) = — d t / du[a(r,u))Zu(r)exp(-iu)t) + c.c. 
27r Jo 

1 Z" 00 

H Vx / dw [/3(r,w)W aJ (r)exp(-ia;t) + c.c] . 

2t io 

The charge and current densities automatically obey the conservation law 

«V + V-j = 0, (27) 

and the remaining two Maxwell equations are also identities because of ([[]). We have 
thus proven that that ([3])-([7j) is the action for macroscopic electromagnetism. 

The electric and magnetic fields are found in the usual manner. Because of (CQ), a 
time derivative of ( )24|) yields the following frequency-domain equation for E(r,w): 

U) 2 

V x [re(r,u;)V x E(r,w)] -e(r, w)E(r, w) = i/x o;j(r,a;), (28) 



j(r, w) = -iwor(r, w)Z w (r) + V x [/3(r, w)W w (r)]. (29) 

The electric field can then be written as 

E(r,i) = — / dw / dV [iwG(r, r', w) • j(r', w) exp(-iwt) + c.c] , (30) 
2vr Jo J 



where the Green bi-tensor G(r, r'u) is the solution of 



uj 2 



V x [ K (r, W )V x G(r, r, w)] - ^-e(r, w)G(r, r, u) = tS(r - r) (31) 

with retarded boundary conditions. Due to ([I]) the magnetic field is given in terms of 
E(r,w) by B(r,w) = -iV x E(r,w)/w. 

3. Quantization 

The canonical formalism [21] is readily applied to the action ()3])-(J7|). A degree of choice 
is present in writing the Lagrangian density £, associated with the action (J3j) by 



5[0,A,X W ,Y W ] = J d 4 x£, 



(32) 



since integrations by parts in ( 132]) change the form of C without affecting the dynamics. 
We take C to be given by the integrands (Ej) (JTJ) , without any integrations by parts. 
The canonical momenta are then 

dC dC f°° 

n * = swT ' nA = ^A) = - £ ° E -/„ *■»»<'■ "J*- < 33 > 
n - = TO) = n -- = flfe = * Y - < 34) 

The familiar complications of canonical QED [2T] are present here: the canonical 
momentum associated with the scalar potential vanishes and the gauge freedom 
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also reduces the number of independent electromagnetic degrees of freedom. These 
difficulties are dealt with |21j by eliminating as a dynamical variable, replacing it 
by its solution obtained from the field equations, and by making a choice of gauge. 
Even after is removed from the list of canonical variables, the vanishing of 11^ is 
responsible [2T] for a constraint on IIa: from the field equation (flOl) and the expression 
for IIa in f|33|) we obtain 

V • n A = 0. (35) 

The general method for quantizing field theories with constraints, such as ( 1351) and 
a choice of gauge, is that of Dirac brackets [21]. For unconstrained theories the 
classical Poisson brackets of the canonical variables are to be replaced by — i/h times 
quantum commutators, but when there are constraints on the canonical variables it is a 
generalization of the canonical Poisson brackets, known as Dirac brackets, that become 
commutators in the quantized theory. The Dirac brackets are computed from the Poisson 
brackets of the quantities that vanish due to the constraints [21]; this requires that the 
constraints be written in terms of the canonical variables. Here ( 1351) is one constraint, 
and a choice of gauge will give another. For free-space QED the Coulomb gauge 

V • A = (36) 

is a convenient choice, but the macroscopic Maxwell equation ( 1231) suggests a 
generalization that has the frequency-domain form 

V- [e(r,w)A(r,w)] = 0. (37) 

In homogeneous media (J37J) reduces to the Coulomb gauge (1361) . but the gauge choice 
(|37|) is not appropriate for inhomogeneous media because in the time domain it cannot 
be written in terms of the dynamical fields and their canonical momenta. As a result, 
Dirac brackets involving the constraint (|37|) cannot be calculated and there is no obvious 
way to proceed with the quantization. We therefore choose the Coulomb gauge (136|) . 
even in the general case of inhomogeneous media. 

After the elimination of as a dynamical field and imposition of the Coulomb 
gauge, the electromagnetic canonical variables are A and IIa, subject to the constraints 
(135]) and (1361) . These are the same constraints as in free-space QED in the Coulomb 
gauge [21] ; an evaluation of the Dirac brackets arising from these constraints shows [21] 
that the quantum canonical operators A and IIa should be defined by the equal-time 
commutation relations 

[A(r,i),n Ai (r',t)] =m Tij (r-r'), 

[A,(r,£), A,(r',t)] = 0, [ft Ai (r,t),n Ai (r',t)] = 0. 

where <5TV/(r — r') is the transverse delta function [2T1 121] 

5 Tij {r - r') = 5 l3 5(v - r') + ( —^—^ " ( 39 ) 
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The canonical operators of the reservoir are subject to no constraints and therefore 
satisfy the commutation relations 



[X ui (r, t), n Xw ,j(r', £)] = ihSij5(u) - u')8(r - r'), 
[X ui {T,t),X u , j (i',t)] = 0, [fl x ^(r, t), H x ^(r', t)\ = 0, 
[Y ui (r,t),fl Yulj (r',t)] = iHSijSiu - u')6(r - r'), 
[Y M (r , t) , Y^- (r', t)] = 0, [ft Yw< (r, t) , fl Yulj (r', £)] = 0. 



(40) 
(41) 
(42) 
(43) 

Both A and IIa are transverse due to the constraints (l36il and ( 1351) and this is consistent 
with the commutation relations fl38|) . From (JT]), f l33]) and fl35|) we see that the transverse 
part of the electric field E is equal to —d t A, while the longitudinal part of E determines 
the scalar potential and is given by the longitudinal part of / °° du a(r, w)X w : 

E = E T + El, (44) 



En 



E, 



-V0 = -e 1 



do; a(r, w)X u 



(45) 



J L 



To write the Hamiltonian we must express the time derivatives of the dynamical 
fields A, X w and in terms of their canonical momenta; this is done for the reservoir 
fields by (I34p . while for the vector potential we obtain from (1331) and (pQ) 



Aa — £qV0 + / da; a(r, w)X 



(46) 



The Hamiltonian density "H is given by 

u = n A -d t A + n x ^ ■ dtX + n Y „ • d t Yu - L (47) 

To express the Hamiltonian H = J d 3 r"H in terms of the canonical variables, all time 
derivatives in (T4T1) must be eliminated using (134"1) and (H6l) . After substitution of ()46l) 
into the first term of the Hamiltonian density (j47j) we encounter a term —IIa ■ V0; upon 
taking the volume integral J d 3 r'H this term becomes, after an integration by parts, 
J d 3 r 0V • IIa, which vanishes by (|35|) . Alternatively, the term — j d 3 IlA • V0 can be 
seen to vanish by noting that V0 is a longitudinal vector, whereas IIa is transverse by 
f l35l) ; the volume integration of the scalar product of a transverse with a longitudinal 
field necessarily vanishes [24] . In the Lagrangian-density term in (T47j) the electric field 
can be expressed in term of the canonical variables by ( 1331) and the magnetic field by 
([!]). The resulting Hamiltonian is 



H 



d 3 r 



IT 



IT 







+ 


-[ 


2 Jo 




1 


+ 








POO 




Jo 



2 
du 



du a(r, w)X a 



+ |(V x A) 2 



ft 2 



+ 1T 2 , 



+ w 2 X 2 +Y,t 



da; a(r, w)X (i 



da/a(r,u/)Xv 



dw/3(r,u;)Y w - (V X A) 



(4£ 
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The Hamiltonian ( |48p gives dynamical equations for the field operators that are 
the quantum versions of the classical equations f fTU|) - f|T3l) . First note that (fit)]) is now 
a statement of the constraint (1351) (see ( 1331) ). rather than a dynamical equation. The 
Hamiltonian ( 148]) and the commutation relations ( 138]) determine the time derivative of 
A as 



d t A = -UA,H] = - 

h e 



n A + ( / dwa(r,w)X 1 , 



(49) 



In this equation the transverse part of du a(r, u)!^ appears because its longitudinal 
part does not contribute the volume integral of IIa • / °° dua(r, w)X w in (]4"8"]) (recall that 
IIa is transverse). The time derivative of IIa is 

9 f n A = ^[II Al if] = - Ko Vx(VxA)+Vx| dw/3(r,u;)Y w . (50) 

Eliminating IIa from (1491) and (I50p gives the dynamical equation for A, which can be 
written in terms of the transverse part of E by ( 145]) : 



-e dtE T -d t ^J dua(r,u)^j = — k V x B + V x J du /3(r, u)Y u . (51) 

This is simply the quantum version of the transverse part of ([H]); the longitudinal part 
of (ITT]) appears here as the time derivative of the second equation in (14"5]K which arises 
from the constraints. The full equation (TIT]) thus holds also in the quantum theory. 
It is straightforward to verify, in a similar manner, that the Hamiltonian (T4H]) gives 
dynamical equations for X w and Y w that are simply the quantum versions of ( [12]) and 
( fl~3]) . respectively. It therefore follows from the results of the last Section that the 
quantum theory specified by the commutation relations ( ]38]) - P3]) and the Hamiltonian 
( 148]) gives the quantum macroscopic Maxwell equations (|23]) -( |24j) for D and H. These 
quantum Maxwell equations will also be explicitly derived in the following Section. 

4. Diagonalization of the Hamiltonian 

We now proceed to diagonalize the Hamiltonian ( T4~8]) . which essentially amounts 
to solving the Hamilton equations for the quantum canonical operators. This 
diagonalization is a familiar procedure from quantization of the microscopic models jH 
COIEJIH], and we follow the method used in [8], where the model was of an inhomogeneous 
dielectric. The field operators in the diagonalized Hamiltonian will be bosonic creation 
and annihilation operators for the eigenmodes of the system. We will show that the 
Hamiltonian can be written 

d 3 r / duhuC{(r,u) -Cx(r,u), (52) 

where C e (r, u) and C m (r, uj) are annihilation operators associated with the electric and 
magnetic responses of the medium, respectively. The diagonalizing operators obey the 
usual commutation relations of bosonic annihilation and creation operators: 

C Xl (r,u),C{ /:j (r',u')] = 6^6(00 - u')6(r - r'), \c M (r, u>), C Vj (r', u')] =0.(53) 



Canonical quantization of macroscopic electromagnetism 11 
From (1521) and (|5"3"|) we obtain the commutators 

C Xi (r,u),H]=hijC x (r,u)), (54) 
which imply that the time-dependent operators are given by 

C A (r, t, oj) = exp(— iu;i)C A (r, lo). (55) 

The canonical operators in the Hamiltonian (1481) must be expressible as a linear 
combination of the diagonalizing operators; for example, it must be possible to write 
the vector potential A and it canonical momentum IIa as 

//»CX> 
dV^ da;[fl(r,r',a;).C A (r',t,a;)+h.c.], (56) 

/poo 
dV/ da;[ff^ A (r,r',a;).C A (r'^,a;)+h.c.], (57) 

where f^r, r', u) and ffi A ( r ) r ' ; are c- number bi-tensors. Relations similar to (023) and 



(157j) must also hold for the remaining canonical operators in (1481) ; thus for and IIx^ 
we require 

X w (r,t) = y; f dV /' 0O da; , [f^(r ) r / ) a;,a; , )-C A (r / ,t,a; , )+h.c.], (58) 

A J J ° 

/rOC 
dV/ da;'[ffi x (r,r',a;,a;0-C A (r',t,a;')+h.c.], (59) 

with analogous formulae for Y w and fly^ in terms of bi-tensorial coefficients fy and fg . 
The bi-tensors and are transverse because of (I36p and (155]) . From (1331) we see 
that fl56|) and (1581) imply an expansion of the electric field E similar to fl56|) . in terms of 
a bi-tensorial coefficient fg related to f£ and by 

f£(r,r» = ^(ry.u)— / du/a(r,u/)f*(r,rV,u;). (60) 

£ o £ o Jo 

To establish the diagonalization we must show that there exist bi-tensorial coefficients 
in f l56|) -f l59|) (and in the corresponding equations for Y w and Tly u ) such that the 
Hamiltonian ( 148]) takes the form (1521) . 

From ( )53|) and ( I56|) - (l59|) it follows that the desired bi-tensorial coefficients must be 
given by commutators of the canonical operators with the diagonalizing operators; we 
give two examples: 

f!(r,r» = [A(r,t),C{(r',t,u)], f*(r, r', u, u') = [X w (r, t), C{(v', t, u% (61) 

The diagonalizing operators must in turn be expressible as linear combinations of 
the canonical operators (the inverse of the relations (|56|) - (l59|) ) and the collection of 
commutators of type (IBTj) . together with (!38l) - P3l) . imply that C A (r,t, lo) has the 
expansion 

C A (r,t,u;)= -i j dV{A(r',t)-f^(r',r,a;)-ri A (r',t)-fl*(r',r,u;) 
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do/ 



X<y (r', t) • f£ (r', r, u/, u) - U x , (r', t) • f£* (r', r, u/, u) 



+%>(r', t) ■ f£(r', r, u/, w) - ft Yw/ (r', t) ■ f Y *(r', r, u>', u)j } 



(62) 



where the asterisks on the bi-tensorial coefficients denote complex conjugation. We 
can now substitute ( 162]) and (14"8"]1 into (151]) and use ([38 ]) -( 143]) to obtain an expansion of 
huC\(r, t, u) in terms of the canonical operators; when the coefficients in this expansion 
are compared with those in the expansion (162]) we find the following conditions on the 
bi-tensorial coefficients: 

POO 

ic^ A (r',r, W ) = «oV x [V x f£(r',r,u;)] - / du'V x [/3(r', u/)f Y (r', r, a;', w)] s (63) 

</ o 



ia;f A (r',r,w) 



to 



CO 



/■oo 

/ du/a(r',u/)f^(r',r,w» 



(64) 



ia;ff) x (r',r,a;',a;) = -a(r' ,u')f* A (r' ,r,cu) + u ,2 f£(r' ,r,cj',cu) 



+—a(r',u') 



du;"a(r',a;")fx(r',r,u;",a;) 



icjfvfr', r, a;', = — f4 (V, r, a;', w 



3 3 3/3 



iwff^r',!-, u/,u;) = -p(j>,u/)V x fX(r',r,w) + u;' 2 f^(r , ,r,w , ,w 
iu;f Y (r', r, a;', w) = 



(r',r,u;» 



(65) 

(66) 
(67) 
(68) 



The last term in (]64"]) is the transverse part in respect of the first index on this bi- 
tensor; the transverse apart appears because this quantity was the coefficient of a 
transverse operator in a volume integral. Note that ( 163]) and (]64"j) are consistent with 
the transverseness of f A and , mentioned after (159]) . We must now solve (l63" ]) - (]68"]) to 
find the required bi-tensorial coefficients. 
Equations (165]) and (166]) imply 

w 2 f£(r', r, to', u) = w' 2 fx(r', r, u', u) 



POO 

f£ A (r',r,u,)+ / dio"a(v'^")^(r\ r ,u",u) 
Jo 



H a(r',u' 



(69) 



where (160]) has been used in the second equality. Equation ( |69]) is essentially the same 
as (I14p . and we write its general solution with the same pole prescription as in (1151) : 

a(r',u') 



to 



a 



(u + iu'0+y 



f^(r', r, u) + 5(u - ^')K(r', r, u 



(70) 



The last term in (170]) is the solution of the homogeneous equation ((J69]) with f£ = 0) 
featuring an arbitrary bi-tensor h x . We do not need to include a term proportional to 
6 (to + to') in (170]) . as is done in ([15]) . because we need only consider positive frequency 
arguments in the bi-tensorial coefficients. From (|67j) and (168]) we find a solution for fy 
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similar to (170]) : 



f Y (r', r, u/, u) 



V x f£(r',r,u) + % - u/)h Y (r',r,u;). 



w' 2 - (w + iw'0+) 2 
Elimination of f£ from ( 163]) and (164]) gives 



13 



(71) 



e w f A (r ,t,uj) 



iuj 



dtu' a(r', u/)f x (r', r, u', u) 



oo 



K V x [V xf^r.w)] 



doo'V x ^(rV^r'^u/,^)] 



(72) 



The integrations in (J72J) can be performed using (I70p and (1711) . in the same way as 
the integrals ffTH]) and (121"]) were evaluated in Section [2j We use the principle-value 
relation (fl~6]) in (170]) and flTTl) . where here the term proportional to 5(u + u/) does not 
contribute as we take only positive frequency arguments in the bi-tensorial coefficients. 
The definitions ([8]) and the Kramers-Kronig relation (Q are also required, and we obtain 

POO 

/ du/ a(r', u/)f A (r', r, u/, w) = a(r', a;)h x (r', r, w) + eofe^, u>) - l]f£(r', r, w), (73) 
./o 

/■oo 

/ da;' f3(r',u)')f$(r',r,u)',uj) = /3(r', w)h Y (r', r, u) - k [k{t' ,uj) - 1]V x f*(r', r, w).(74) 
./o 

Equations fl60|) and ([Mjl show that 



k^(r / ,r, W ) = f^(r',r,u;) + - 



UQ 



J L 



which implies 

icoV x f A (r',r,u;) = V x f£(r',r,u). 
Insertion of (1751) . (174]) and (175]) in ([72]) and application of (176]) produces 



(75) 



(76) 



V x Kr»V x f£(r',r,a;)] - ^(i»f A (r', r, u; 



/i w a(r', cu)h x (r', r, u) + i/i uV x [/3(r', u;)h Y (r', r, a;)] 



(77) 

This equation is essentially the same as the wave equation (I2"8"j) - (l2"9"j) for the electric 
field in the frequency domain, and its solution can be written as in (130]) using the Green 
bi-tensor (13~T1): 



f A (r',r, W )=/i J dV'G(r',r",a;)-s A (r",r,u;), (78) 

s A (r', r, to) = co 2 a(v', w)h x (r', r, to) + iwV x [/3(r', w)h Y (r', r, w)]. (79) 

All of the bi-tensorial coefficients can now be written in terms of h x and h Y ; to 
avoid overly lengthy expressions we will write all the coefficients in terms of h x , h Y and 
f^, bearing in mind that fg is itself determined by h x and h Y through (1781) - (1791 . The 
right-hand sides of (163]) -( 168]) are written in terms of h x , h Y and fg using ( 170]) . ( 17X1) . 
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(J73J, dZU), (En}, dZZD> with the results 

-£ £(r', w)f^(r', r, u) - a(r', u;)h£(r', r, uj), 



fn A (r',r, W 



fE(r , ,r,w)] T , 



UJ 



(80) 
(81) 



f n x ( r '' r ' a;/ ' a; ) = - a ( r 'iu')f B (r',r,u) - iwh x (r', r, w)5(u; - u/) 



2co 



1 



+ 



1 



u/ — uj — iuO + u/ + uj — iwO 



fE(r',r, W ) 



i 



^(r',r,uj',uj) = -^ x (r',r,uj',uj), 

UJ 



f* (r', r, uj', u) = -^/?(r', u/) V X f*(r', r, uj) - iwh Y (r', r, uj)8{uj - uj') 



(82) 
(83) 



A /„/ 



u'f3(r',uj') 
2u7 2 



1 



1 



uj' — to — iuO + uj' + to — iwO + 



f*(r , ,r, £ y,o;) = -f5 Y (r',r,a; , ,w 



J4) 



(85) 



Equation ( 1771) shows that the right-hand side of ( 1501 is transverse, a fact that is used 
in deriving (180|) and ( 15Tj) . The pole prescription in ( 170"]) and (171 p has been rewritten in 
( 1521 and ( 1541 as shown in the first equality in (ITBl . 



All that remains is to deduce the bi-tensors h x and h Y ; once these have been 
specified the bi-tensorial coefficients are found from ( 1751) -( 1851) . The conditions that 
determine the form of h x and h Y follow from substituting the expansion ( 162]) and its 
hermitian conjugate into the commutation relations ( 153]) . When this substitution is 
made in the first of (1531) . and the relations (1501 . (15T1) . ( 1531 and ( 1551) are used, we obtain 
the requirement 

dV {--, [s e*(r", uj)f^ ki (v", r, uj) + a{v" , c)^(r", r, uj)] f£ kj (r", r', u/) 



UJ 

i 



/&(r", r, uj) e e(v", uj')f£ kj (r", r', u/) + a(r", uj')h* kj (v", r', a;') 



i i 

- + - 

uj' uj 



du" 



f^ kt (r",r,uj",uj)fl kj (r",r',uj",uj' 



i I II H \ e\ I II I II I 



—ihSij5\\'5(uj — uj')5(t — r'). 



(86) 



The quantities in square brackets in the first two lines of (1561 are transverse because of 
(177]) ; they therefore project out the transverse parts of the bi-tensors and fg* with 
which they are combined in a scalar product. For this reason the transverse restriction 
on fg' and f^*, which was present because of (1811) . has been removed in (156"]) . To find 
the condition on h x and h Y that results from (1561) . we must insert (1521) and (1541 and 
perform a simplification that makes use of (1771) ; the details of this tedious calculation 
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are described in Appendix A and the result is 



d d r 



h$ ki (r", r, u)h^ kj (v», r', u ) + K* ki (r", r, u)f4 kj (r", r', w ) 
ft 



2lu 



SijS X yS(r - r') 



57) 



We must also impose the condition arising from the substitution of (162]) into the second 
commutation relation in (I55|) . But this condition is found to be automatically satisfied 
when fl80l - fl85l) and ([77]) are utilized, and we relegate the details of this calculation 



to Appendix B 



The diagonalization ( 152]) is achieved once we find bi-tensors and h Y that satisfy 
( ]87j) . The diagonalization operators are of course not unique since, for example, a unitary 
transformation C^(r, oS) — > UijC\j(r,uj), Uj k Ukj = does not change the form ( 152]) . 
and similarly a unitary transformation can be applied to and h Y without affecting 
the condition (1571) . A suitable solution of (1ST]) is 

V2 / \ 1/2 



^i(r',r,a;) 



5Ac ^' <5(r " r/) ' ^ (r '' r ' W) = (2^) <"^ r - r ')' ( 88 ) 

which completes the diagonalization procedure. 

We can now write the electromagnetic field operators in terms of the diagonalizing 
operators C\(r,t,u>). The bi-tensor s A in (!78 l) -( !79l) is seen from ([88]) to be 



s c (r', r, oj) = uj 1 



he , 
ei(r,u) 

71 



1/2 



s m (r',r,w) = iuV x 



16 (t - r) 

1/2 



7T 



/ / X 

«i(r ,uj) 



15(r - r') 



(89) 



(90) 



which determines through ( 178]) the bi-tensor fg. Having obtained f^, the electric-field 
version of (1561) together with (153]) show that 



Z7T 



do; / dV 



J(r,wj 



-27riu; 



iwG(r, r', u) ■ j(r', w) exp(— icut) + h.c. 



(91) 







nl/2 



7T 



-£i(r,o;) 



C e (r,w) + 2ttV x 



h,Kr 



71 



Ki(r,u)) 



nl/2 



C m (r,wU.(92) 



The current operator defined in ( 192]) is easily seen from ([53]) to obey the commutation 
relations 



j(r, W ),jt(rV) 



j(r,w),j(r',u/) 



47rM(a; - u') I u 2 e ei(r, u)15(r - r') 



+ k V x 



0. 



^-ki(v,ui) tS(r - r')y/-Ki(r 



' u') 



x V' 



(93) 
(94) 



where the notation x V denotes a curl with respect to the right-hand index, so that 
V(r) x V= V x V(r) for a vector V(r). The transverse and longitudinal parts of E(r, t) 
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are related to the electromagnetic potentials by ( 145)) and the magnetic field operator 
can be found from 

B(r, u) = -iV x E(r, u) /u, (95) 

which is an identity in terms of the electromagnetic potentials. 

Equations (191"]) and (I92j) are the quantum versions of the classical results (130]) and 
( 1291) . The diagonalizing operators C e (r, u) and C m (r, uj) are thus seen to be closely 
related to the classical fields Z^(r) and W u (r) in Section [2j It is easy to show that (1911 
is the solution of the quantum macroscopic Maxwell equations. From ( 1911 and ( 1951) . 
and the defining equation (|31~!) of the Green bi-tensor, it follows that 

VxH-^D=j, (96) 
with use of the definitions (pED and ([22]). Similarly, (|9"T1). (Tl9]) and (151]) imply 

V • D = a, (97) 
where in the frequency domain the charge density is 



a(r, u) = — V ■ j(r, u) = -2ttV 

100 



he 



71 



1/2 



C e (T,u)\, (98) 



so that the quantum version of the conservation law (!2Tj) holds. 

In the phenomenological approach to macroscopic QED [5] [6], a quantum noise 
current j(r, uj) obeying the commutation relations (|93j)- (l94j) is postulated, and this 
current together with its associated charge density a(r, uj) = V • j(r, oj)/(ioj) are 
inserted into the macroscopic Maxwell equations ( l96]) -( l9~Tj) . This leads to quantum 
electromagnetic field operators that are then shown to satisfy appropriate commutation 
relations as a consequence of (l93]) - (|94l) . The noise current is further decomposed 
into bosonic operators C e (r, a;) and C m (r,o;), as in ( |92]) . and ( 152]) is then the 
required Hamiltonian for these bosonic operators. Here we have canonically quantized 
macroscopic electromagnetism, proved that C e (r,co>) and C m (r, uj) in (|9T ]) -( l92]) are 
the diagonalizing operators of the Hamiltonian, and derived the current-density 
commutation relations (|93])- (|94l) . 

An interesting issue raised by Huttner and Barnett [1] is the possibility of negative- 
energy eigenstates in their microscopic model if the coupling constants are too large, 
since this can happen in the coupling of harmonic oscillators. Huttner and Barnett [1] 
found only positive frequencies in the diagonalized Hamiltonian, without restrictions 
on the size of the coupling constants in their theory, in which only the transverse part 
of the vector potential is quantized. This conclusion was supported by verifying that 
operator expansions analogous to (|56 ]) -( |59]) can be inverted to give expansions analogous 
to ( 162]) . with only positive frequencies present in the expansions. Here we have similarly 
found that only positive frequencies are required in the diagonalized Hamiltonian (152]) . 
The relation of these results to the case of general coupled harmonic oscillators and the 
occurrence of negative energies may be interesting to investigate. 
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5. Conclusions 



The absorption of electromagnetic energy in media is a necessary consequence of the 
restriction to retarded solutions of Maxwell's equations, which leads to the Kramers- 
Kronig relations. This means that the canonical quantization of electromagnetism in 
media cannot proceed by treating only the electromagnetic fields as dynamical variables, 
unless restrictions are made in the frequency range of applicability. By including a 
reservoir in addition to the electromagnetic fields, we developed macroscopic QED using 
the standard canonical formalism of quantum field theory. The theory applies to any 
linear, inhomogeneous, magnetodielectric medium with dielectric functions satisfying 
the Kramers-Kronig relations. The starting point of the phenomenological approach to 
macroscopic QED was derived from the underlying canonical theory. Important results 
such as the Lifshitz theory [I7J [TH1 EE] of vacuum and thermal electromagnetic forces 
can therefore be provided with a canonical foundation. 
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Appendix A. 

Here we describe how (!86|) is simplified to (187|) . The first step is to insert ( )82|) and ([8j 
into the terms in (1861) that depend on and f^ Y and their complex conjugates. We 
must note that all frequency arguments in f l86|) are positive, since they are all positive 
in f )53|) . which leads to fl86|) . Consider first the term 

/ dV {- (1 + 1) ^d^/^^.r.t/.^^y.^ty)^ (A.l) 

in ( 1861) . Insertion of ( 1821) into (1 A. II) leads to products of the pole prescriptions in 
(J52|) with their complex conjugates; these various products can be rewritten using the 
following identities: 

1 < 1 * '(A,) 

(A.3) 
(A.4) 

(lo" + lo + iO+)(w" + to' - i0+) ~ u - lo' + 2i0+ V^" + fa/ - i0+ to" + lo + i0+ J 

By expanding the quantities in brackets in these formulae in terms of principle values 

and delta functions, as in ( Tl6|) . and similarly expanding the other pole prescriptions in 



(w" -oj + iO+)(u/' - a/ - i0+) uj-uj' - 2i0+ V^" - w + i0+ fa/' - u' - iO" 
1 1/1 1 



(lo" -uj + i0+) (oo" + u' - i0+) w + w' - 2i0+ V w" - w + i0+ u" + lo' - i0~ 
1 1/1 1 



{lo" + lo + iO+)(fa/' - fa;' - i0+) lo + lo' + 2i0+ V^" - w' - i0+ lo" + lo + iO H 
1 1/1 1 
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(lA.ip that involve oj", the integration over oj" can be evaluated. As all of the frequencies 
u, {J and oj" are positive, any delta function whose argument is a sum of two frequencies 
does not contribute. Making use of the Kramers- Kronig relation in performing the 
u" integration, and after straightforward but lengthy manipulations, we obtain for (lA.lj) 



-i / d d r 



3„'/ 



— + m5(uj — oj') 



oj — oj 

xf^(r",r,oj)f^(r",r',oj 
1 



e [e*(r",u))-e(r",u)')) 



+ 



u) — OJ 
1 



- + m5(uj — oj') 



- + m5{uj — oj') 



a{v'\u)')f^Uv\v^)hi kj {v'',v\u)') 
a(r",u))h$ ki (r^T,u)fg kj (r",r',u)') 



OJ — OJ 

+2oj5(oj - u)')h£ ki (r", r, u))h£ kj (v", r', u 
Proceeding similarly with the ffj -dependent term in fl86l) yields 



(A.6) 



3„" 



d d r 



-i / d 3 r" 



l l 

- + - 

U)' OJ 



doj''f^ kt (r'',r,oj'\oj)fl kl (r'',r',oj'',oj') 



WW 



P + m5(oj - oj') 



OJ — OJ 



K [k(t",u)') - K*(r",u)] 



x(V" x ^) ki (v",v,u))(V" x i£) kj (r",r',oj" 



i 

U) 

i 

"u7 



OJ — OJ 



- + i7c5(oj — oj') 



/3(r",u/)(V" x f**) ki (r",r,u))h* kj (r'',r',u)'" 
0(r",u))h^ ki (v",v,u))(V' x f£) kj (r",r',u)') 



P h md(oj - oj') 

oj — oj' 



+2oj5{oj - w')4(r", r, cu)h* kj (r", r', oj') } . 



(A.7) 



The results ( 1A.6I) and ( 1A.7j) must now be used in ( 1861) . Combining the terms in the first 
two lines of ( |86|) with ( 1A.6|) produces the quantity 



i / dV <! 5 



OJ 



oj'(oj — a;') 



e*(r",u))-P 



OJ 



oj(oj — a;') 



— e(r , u/) 



+i7r[6*(r", w) - e(r", a/)]5(u; - a/) 



/E^(r",r,u;)/ E - (r",r',w 



A' / // / , ,/\ 



EfcjV 



a;' 



a; (a; — u/) 



i7r5(w — a;') 



a(r",a;0/ E ^(r",r,w)^ 7 (r",r',w 



+ 



OJ 



— + i7r5(w - oj') 



oj'ioj — oj') 

+2u)S(u - oj')h^ kt (v", r, w)Z#y(r* r', u)')\. 



a(r",w)^(r",r,w)/ E A l ? (r",r',u;') 



(A.J 
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We must now add flA.8[) to ( )A.7j) and thereby obtain the left-hand side of (15BT) . Before 
performing this addition we make an integration by parts in each of the terms in (I A. 71) 
containing a curl. It is then straightforward to see that because of ( 177|) most terms in 
the sum of (1A.8j) and (1A.7I) cancel. The only terms that survive are the last term in 
(TAH and the last term in flATfl) . We thus obtain (1571). 



Appendix B. 

Here we describe the result of substituting fl62|) into the second commutation relation 
in ( )53l) . This calculation is of course very similar to the substitution of ( l62l) into the 
first commutation relation in ( 153|) . which simplified to the condition ( 1H7|) . In this case, 
instead of the condition ( 186|) on the bi-tensorial coefficients, we obtain 

dV \-- f [e e*(v" : u)f*Ur", r, oj) + a(r", oj)h$ ki (v", r, oj)] /^*(r", r', oj') 



oj 



+ r, oj) e e*(r'>, a;')/^(r", r', u/) + a(r", (r", r', oj') 



i i 

u/ UJ 



du" 



f^ kt (r'',r,uj'',uj)f^ kj (v'',r',uj'',uj 



A'* /__// „/ , ,// , ,/\ 



+ /£*(r", r, u/', c)/^.(r", r', w ", a/) = 0. (B.l 



Proceeding as in Appendix A , we substitute (1521) and flM|) into fIB.ip . Consider first the 
term 

i i 



3„// 



duj''f^ kl (r'',r,uj",uj)f^ kj (r"y,uj",uj'^ 



(B.2) 



(w" - oj + i0+)(u" -uj' + i0+) w - w' W - w + i0+ w" - u/ + i0+ 
1 1/1 1 



Insertion of (1821) into (1B.2I) leads to products of the pole prescriptions that can be 
rewritten using the following identities: 

1 1/1 1 

(B.3) 
(B.4) 
(B.5) 



(uj" -uj + i0+)(uj" + uj' + i0+) uj + uj' \uj" - uj + i0+ oj" + uj' + i0+ 
1 1/1 1 



(uj" + uj + i0+)(uj" -uj' + i0+) uj + uj' \uj" - uj' + i0+ uj" + uj + i0+ 

\ = _J_ ( 1 (B 6) 

(uj" + uj + i0 + )(uj" + uj' + i0+) uj-uj' \uj" + uj' + i0+ uj" + uj + i0+ / 1 ' J 

In the terms l/(u — uj') on the right-hand sides of identities ( IB. 31) and (1B.6I) there is no 
prescription for dealing with the pole at uj = oj'; this is not a problem however, because 
when ( 1B.3[) -( jB~6l) are used in ( IB. 21) these factors of 1/(uj — oj') get combined with other 
terms in such a manner that they do not give rise to a pole. Bearing in mind that all 
frequencies are positive we obtain for (IB. 21) 



3„" 



d d r 
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-i / dV 



- f e [e V, u>) - e*(r", uj')) / E ^(r", r, u)f^(f, r>, uj') 



+ 



UJ + UJ 

1 



A'* /„// / , ,/\ 



UJ + UJ 

1 



uj + uj 



- i a(v'',u)h^ i (v'',v,u)f^ j (v'',v' 



UJ 



(B.7) 



Note that there is no term quadratic in h x in (jB.7j) . in contrast to the analogous result 
( 1A.6I) . This is because in the case of ( IB. 71) the term quadratic in appeared with a 
factor (u — U)')5(ou — uj') = 0. The fn Y -dependent term in (IB. 1[) similarly yields 



-i / d d t 



3„" 



ujuj\u + u/J 

x(V" x ^*) ki (r",r,u)(W x f£'%(r", r', u/) 



— - — -/3(r",u/)(V" x f£%(r* r, u,)/^", r', uj') 

UJ\UJ + OJ J 

- 1 /j(r",c)^(r-,r,a;)(V- x f£'%(r", r', uj') 
a;' (a; + u/J 



(B.8) 



Combining the terms in the first two lines of (IB.lj) with ( 1B.7|) produces the quantity 



- i / dV <! e 



UJ 



uj'{uj + u/) 



u;(u; + a;') 



A'* /„// „/ , ,/\ 



UJ' 



A'* /„// _/ , ,/\ 



U){uj + UJ 1 ) 



+ 



— a(r", w)/i^(r", r, w)/^*- (r", r, a/} 



(B.9) 



+ uj') 

We now add (1B.9|) to (IB. 81) to obtain the left-hand side of (IB.lj) . An integration by parts 
in each of the terms in (1B.8j) containing a curl and use of (ITT)) shows that the result of 
this addition is zero. The condition ( IB.lj) therefore gives no restriction on the bi-tensors 
hx and h Y . 
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